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[i] We present in situ and ground-based measurements of the ratio k of ion cyclotron 
angular frequency to ion-neutral momentum transfer collision frequency to investigate 
ionosphere-thermosphere (IT) coupling in the auroral E region. In situ observations were 
obtained by NASA sounding rocket 36.234, which was launched into the nightside 
E region ionosphere at 1229 UT on 19 January 2007 from Poker Flat, AK. The payload 
carried instrumentation to determine ion drift angle and electric field vectors. Neutral winds 
were measured by triangulating a chemical tracer released from rocket 41.064 launched 
two minutes later, k is calculated from the rotation of the ion drift angle relative to the 
E-cross-B drift direction in a frame co-rotating with the payload. Between the altitudes of 
118 km and 130 km k increases exponentially with a scale height of 9.3 ± 0.7 km, 
deviating from an exponential above 130 km. k = 1 at an altitude z 0 of 1 19.9 ± 0.5 km. 

The ratio was also estimated from Poker Flat Incoherent Scatter Radar (PFISR) 
measurements using the rotation of ion velocity with altitude. Exponential fits to the PFISR 
measurements made during the flight of 41.064 yield z 0 = 115. 9 ± 1.2 km and a scale 
height of 9.1 ± 1.0 km. Differences between in situ and ground-based measurements show 
that the E region atmospheric densities were structured vertically and/or horizontally on 
scales of 1 km to 1 0 km. There were no signs of ionospheric structure in ion density or ion 
temperature below scales of 1 km. The observations demonstrate the accuracy with 
which the in situ and PFISR data may be used as probes of IT coupling. 
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1. Introduction 

[ 2 ] We present results from a new technique for measuring 
ionosphere-thermosphere (IT) momentum coupling in the 
E region ionosphere from a spinning spacecraft, and we 
compare these results with simultaneous observations from 
an incoherent scatter radar (ISR). Collisional coupling 
between the ionosphere and neutral atmosphere represents an 
important pathway for the transport of mass, momentum, and 
energy in the space environment. Collisions determine where 
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field-aligned currents in the high latitude ionosphere close and 
thus where magnetospheric energy is deposited [e.g., Thayer, 
2000]. Numerous satellite observations have indirectly dem- 
onstrated that the neutral atmosphere responds to ionospheric 
forcing at high latitudes, especially in the auroral and magnetic 
cusp regions [Emery et al, 1985 \Heelis, 1987; Knowles et al., 
2001; Liihr et al., 2004; Clemmons et al., 2008; Prolss, 2008; 
Sakanoi et al., 2009; Ritter et al., 2010], The resulting varia- 
tions in neutral density affect low-perigee satellite orbits and 
consequently impact the accuracy of satellite orbit prediction 
and tracking operations [e.g., Storz et al., 2005]. 

[ 3 ] Incoherent Scatter Radar (ISR) can probe IT coupling 
processes [Dougherty and Farley, 1963], and many ion- 
neutral collision frequency estimates have been reported 
using the ISR technique [e.g., Reese et al., 1991; Nygren, 
1996; Davies et al., 1997, and references therein]. Equally 
significant is that ion-neutral collisions play a crucial role 
in deriving neutral winds, electric current densities, Joule 
heating rates, and neutral atmospheric properties from ISR 
measurements [Alcayde et al, 1972], In some cases the ion- 
neutral collision frequency is derived directly from fitting the 
ISR power spectrum or autocorrelation function to a model, 
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Figure 1 . Illustration of 36.234 experiment payload (not 
to scale) showing the Suprathermal Ion Imager (SII) and its 
coordinate system, the forward Ionization Gauge (IG), and 
the Electric Field Experiment (EFE) booms. (The payload 
bus and EFE booms have been foreshortened for the purpose 
of illustration.) The payload symmetry axis was aligned to 
within a few degrees of the geomagnetic field vector (Bo) 
by an attitude control system. The payload spin period was 
1.6 s. 

which typically involves assuming the ion and electron 
temperatures are equal and/or equal to the neutral tempera- 
ture. An alternative technique uses measurements of the 
altitude variation of ion drifts [e.g., Nygren et al., 1987, 1989; 
Davies et al., 1997]: Ions transition from drifting with the 
neutrals at low altitudes (collisional regime) to drifting at the 
E x B velocity at high altitudes (collisionless regime). As 
the neutral winds usually are not measured, this method 
typically assumes they are negligible. Uncertainty about the 
validity of the assumptions that go into each technique makes 
automation difficult. So, despite substantial progress in this 
area, collision frequencies usually are not quoted in ISR 
observations. Combining ISR measurements with simulta- 
neous in situ observations provides a means of examining 
and validating those assumptions, and this is one element of 
our study. 

[ 4 ] In contrast to the many ground-based observations 
cited above, only two direct in situ observations of ion- 
neutral coupling in the E region have been reported in the 
literature. The pioneering study by Watanabe et al. [1991] 
used sounding rocket observations to demonstrate ion 
demagnetization in the E region ionosphere, although neutral 
winds were not measured. Progress was made recently by 
Sangalli et al. [2009] who used ion drift, electric field, and 
neutral wind measurements obtained by a pair of sounding 
rockets from the Joule-II campaign to self-consistently esti- 
mate the height z 0 of the boundary between collisional and 
collisionless ion drift to be 1 18.0 ± 0.3 km during an auroral 
breakup near Poker Flat, AK, on 19 January 2007. 

[5] In this paper we analyze ion drift, electric field, and 
neutral wind measurements obtained by a second pair of 


rockets (36.234 and 41.064) from the Joule-II mission, 
during the same auroral breakup, and derive a height profile 
for the ratio k of ion cyclotron frequency to ion-neutral 
momentum transfer collision frequency. We validate the in 
situ results by comparing scale heights in k with neutral 
atmosphere pressure scale heights measured by an ionization 
gauge. We then compare these measurements, and those of 
Sangalli et al. [2009], with measurements from the Poker 
Flat Incoherent Scatter Radar (PFISR) to assess the accuracy 
of the radar observations of IT coupling. The instruments and 
observations are described in section 2. Section 3 presents the 
techniques and results, which are discussed and summarized 
in section 4. Appendix A provides details on both the error 
analysis and the ion-neutral momentum transfer collision 
frequency model calculations. 

2. Instrumentation and Data 

[6] NASA sounding rocket 36.234 was part of the Joule II 
mission to study frictional heating in the E region ionosphere 
[Sangalli et al., 2009; Hvsell et al., 2008], The rocket laun- 
ched from the Poker Flat Research Range at 1229 UT on 
19 January 2007 during a period of relative calm in both the 
aurora (quiet, stable arc) and geomagnetic activity (K p ~ 2 + 
to 3—, a p ~ 10 nT, and DST ~ —10 to —15 nT in the hours 
surrounding the launch). An attitude control system (ACS) 
was used to orient the payload rotational axis to within a few 
degrees of the local geomagnetic field direction and to 
maintain the payload spin period at 1.6 s. The payload car- 
ried multiple instruments to an altitude of 190.2 km, 
including a Suprathermal Ion Imager (SII), which takes 
snapshots of the ion energy/arrival angle probability distri- 
bution function at a rate of 125 distributions per second. The 
SII was mounted at the end of a 1 m boom on the rocket 
payload as illustrated in Figure 1. Details related to the SII 
principle of operation are given by Knudsen et al. [2003]. Of 
relevance to the present study is that the SII analyzer projects 
an image of the differential directional energy flux of ions 
arriving within ±5° of the plane of its cylindrical entrance 
aperture onto a microchannel plate intensifier (Figure 2). 
Because the payload spins about its symmetry axis (see 
Figure 1), the total detector counts are modulated as the 
entrance aperture sweeps across the ion ram direction. There 
is a peak in number of ions entering the aperture when the 
ram vector lies within the aperture field of view. On this 
flight there was a fracture in the MCP detector which pre- 
cluded obtaining estimates of the ion velocity vector. As will 
be shown in section 3, we estimate the direction of the ion 
flow vector by determining the relative spin phase at which 
the integrated ion flux is a maximum. See Drakou et al. 
[1997] and Knudsen and Wahlund [1998] for other appli- 
cations that use the total flux from a spinning probe to 
determine plasma properties. Figure 3 shows an ion distri- 
bution image that is representative of the data obtained on 
this flight. It is evident from watching a movie of the image 
time series that the MCP intensifier is cracked, and this 
distorts the distribution of signal within the image. Therefore 
we have selected for analysis the lower-right hand side of 
each image, where the MCP is operating nominally. The 
selected data correspond to altitudes below 155 km on the 
descent stage of the rocket’s trajectory. The circular fiducials 
represent kinetic energies of 2 eV and 5 eV, assuming singly 
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Figure 2. Illustration of the Suprathermal Ion Imager 
principle of operation. Figure 2 represents a 2D section 
in the x-z plane of the coordinate system in Figure 1. After 
Burchill et al. [2010]. 


ionized ions. A bias of —1.9 V relative to payload ground 
was applied to the sensor skin, and the peak signal lies just 
outside the inner fiducial, consistent with a slightly negative 
payload floating potential. The integration period for the ion 
distribution images is 8 ms. 

[ 7 ] Rocket 36.234 also carried a multiple-sensor ionization 
gauge in order to measure the state of the ambient thermo- 
spheric gas. This instrumentation consists of five pressure 
sensors oriented such that the ambient flow is sampled at five 
angles of attack. The sensors employ miniature Bayard- 
Alpert ionization gauges [ Bayard and Alpert, 1950] and 
modem, high dynamic-range electrometers for sensitive 
pressure measurements. Each gauge measures the pressure 
inside its own accommodation chamber. One sensor, the 
forward sensor, was located on the forward payload deck and 
oriented with its entrance aperture normal aligned with the 
spin axis. The purpose of this sensor was to enable temper- 
ature and density profiles to be determined via the standard 
“hypsometric” method [e.g., Horowitz and LaGow, 1957]. 
For this study, we use the neutral pressure measurements on 
the upleg to validate the re scale height measurements. 

[8] The DC and AC vector electric fields were measured 
using spherical sensors with embedded pre-amps extended 


on 5.2 m tip-to-tip articulated fiberglass booms that unfolded 
after the payload had attained sufficient altitude in the lower 
ionosphere. Since the payload spin axis was oriented along 
the magnetic field direction, the two-dimensional electric 
field measurement completely parameterized the electric 
fields perpendicular to the magnetic field direction. The DC 
electric field data were sampled with 1 8 bit A/D converters 
with a time resolution of 1 msec. The v pa yioad x B fields were 
subtracted from the measurements and the resulting DC 
electric fields were then rotated into geomagnetic coordinates 
using the attitude solutions provided by the on-board gyro- 
scope. Electric field data gathered during this flight are ana- 
lyzed and discussed by Hysell et al. [2008]. Data used for the 
present analysis were obtained from fits of one component at 
twice per spin (1.25 Hz) to reduce measurement noise. 

[ 9 ] A second rocket (4 1.064) launched at 1231 UT released 
a luminescent tracer (trimethyl aluminate) which was trian- 
gulated optically by ground stations to determine the zonal 
and meridional components of the neutral bulk velocity 
vector. In the following analysis we use data from the upleg 
of 41.064, as measurements were not available from the 
downleg. 

[ 10 ] The in situ measurements were complemented by 
ionospheric observations from the Poker Flat Incoherent 
Scatter Radar (PFISR) facility. PFISR operated in a mode 
with seven beams, as illustrated in Figure 4 and described 
by Heinselman and Nicolls [2008] and Hysell et al. [2008, 
2009]. PFISR transmissions consisted of interleaved long 
pulses and alternating codes. The long pulse and alternating 
code data were integrated for two and five minutes, respec- 
tively, prior to performing incoherent spectral analysis. The 
multibeam, pulse-to-pulse steering capability of PFISR 
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Figure 3. Representative ion image at T + 294.36 s. The 
dashed circular fiducials represent kinetic energies 2 eV 
and 5 eV, assuming singly ionized particles. Pixel count rate 
is represented linearly by the gray scale. The signal peak lies 
slightly outside of the 2 eV circle. This is consistent with 
negative payload charging of a fraction of a volt relative to 
the plasma. 
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Figure 4. PF1SR beam geometry and rocket trajectories in 
altitude-adjusted corrected geomagnetic coordinates, (bot- 
tom) Magnetic longitude vs. magnetic latitude and (top) 
altitude vs. magnetic latitude. The radar beams (dotted lines) 
are numbered 1 through 7, and beam 7 points anti-parallel to 
the geomagnetic field line. Curves denote sounding rockets 
36.234 (dashed) and 21.138 (solid). 


allows vector drifts, and corresponding inferred electric 
fields and neutral motions, to be derived. 


3. Techniques and Results 

[n] Under the steady state guiding center approximation, 
ions drift with the average velocity v,x in the direction per- 
pendicular to the geomagnetic field vector according to 


^ + S± x *) + x h/B ' 


where re = is the ratio of the average ion cyclotron 

angular frequency to the average ion-neutral momentum 
transfer collision frequency, E is the electric field vector, B 
and b are the magnitude and direction of the geomagnetic 
field, and u is the average drift velocity of the neutral 


atmosphere. Equation (1) is a special case of equation 
(5.103) of Schunk and Nagy [2009], for which ion pressure 
gradients and gravity are negligible. All neutral species are 
assumed to drift with the same velocity, as are all ion spe- 
cies. In what follows we describe how equation (1) is used 
to determine re from the in situ measurements, and describe a 
method to estimate the effect of collisions using 1SR 
observations. 

3.1. In Situ Method and Results 

3.1.1. Ion Ram Demagnetization 

[ 12 ] Sangalli et al. [2009] calculated re from the ratio 
| Vo. — u±\/\E ± /B — v,_l x b | (see equation (1)) using mea- 
surements of the electric and magnetic fields, the ion drift, 
and the neutral wind from rockets 21.138 and 41.065. In the 
absence of vector ion drift measurements on 36.234, we 
solve equation (1) for Vjx to calculate the deflection 9 of the 
ion ram vector due to ion-neutral collisions assuming a 
particular value of re, and find that re for which the calculated 
deflection matches the observed deflection. Figure 5 illus- 
trates the geometry of the calculation. From measurements 
of the payload velocity (v p ±), the neutral velocity («x), and 
the electric field (E±), we calculate the ion velocity for each 
of a set of re values running from io -L500 to 10 L ~ 00 in fac- 
tors of IO 0 001 . From these we calculate the expected angle 9 


8 = arccos 


( U'X ' ^ExB \ 

VluxlliWly 


(2) 


between the observed ion ram vector v- ± = v,x — v p i_ (ram) 
and the rammed E-cross-B drift vector v' ExB = v Ex b ~ Vp± in 
the frame of reference co-rotating with the spinning payload. 
As seen in the payload frame, nonzero values of 9 represent 
a “demagnetization” of the ion ram vector due to collisions 
with neutrals. The true re is obtained where 9 matches the 
observed deflection angle, 9 m . We perform the calculations 
in a plane perpendicular to the magnetic field vector calcu- 
lated from the eleventh generation International Geomag- 
netic Reference Field model [Finlay et al, 2010]. 

[ 13 ] Figure 6a shows a representative time series of the S1I 
count rate obtained from a sum over the right half of the SII 
images. Figures 6b-6d show selected spin profiles with four- 
parameter Gaussian fits (dashed curves) at three different 
altitudes of the rocket. The gray curves in each panel are 
components of the v ExB projected onto the SII entrance 
aperture. The time difference At between the peak in the SII 
signal and the peak in the projection of v ExB is converted to a 
ram deflection angle (the angle between the ion drift vector 
and the E x B vector in the payload frame of reference) 
according to 

8 m = (A t/T) x 360° (3) 

where T= 1.6 s is the spin period of the payload. We use the 
fact that 9 m = 0° above 140 km to calibrate the angle 
between the SII boom and the electric field x-axis to within 
0.25°. Figure 6a is the count rate profile at an altitude of 
~149 km, where the ions are E-cross-B drifting. Figure 6c 
shows the profile measured at ~125 km, showing a clear 
demagnetization of the ion ram vector. The E-cross-B drift 
can either lead (as shown) or lag the SII signal, depending 
on the relationship between E ±_ , v p ± , and u± . Figure 6c 
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Figure 5. Illustration of deflection of ion ram vector due 
to ion-neutral collisions (“ion ram demagnetization”). 
From measurements of the payload velocity (v p ), the neutral 
velocity (u±), and the electric field (E±), we calculate the 
expected angle 9 between the observed ram vector and 
the E-cross-B vector v ' ExB for a range of re values. We obtain 
the true re where 9 matches the observed deflection angle, 

^measured* 

again shows significant demagnetization, at an altitude of 
~ 1 1 0 km. 

[ 14 ] Figure 7 shows the altitude profile of the measured 
deflection angles. Between 140 km and 150 km the standard 
deviation of the deflection angle is 0.25°, and we use this as 
an estimate of the uncertainty in the deflection angle mea- 
surement, A 9 m . 

[is] We use the calibrated ram deflection measurements 
between 118 km and 140 km to determine re according to the 
technique described above. Analysis of the uncertainties in re 
are presented in section Al. The angular deflection is 
anomalously large for the observed neutral wind and electric 
field at altitudes below 117 km, and we are unable to obtain 
measurements of re in that region. Figure 8 shows the 
resulting height profiles of re and v in between 118 km and 
140 km. To identify the height z 0 at which the ion-neutral 
momentum transfer collision frequency matches the ion 
cyclotron angular frequency, we fit a straight line to the 
model 


hire = (z-z 0 )/H k , (4) 

over the height range 117-131 km where the profile is 
approximately exponential. The result is the straight line in 
Figure 8a, where z 0 = 1 19.9 ± 0.5 km and the corresponding 
re scale height is H K = 9.3 ± 0.7 km. For comparison we 
have calculated the theoretical ion-neutral collision fre- 
quency (dashed curve) for model neutral densities from 
NRL-MSISE-00 [Picone et al, 2002] and IRI 2007 \Bilitza 
and Reinisch, 2008] using non-resonant polarization attrac- 
tion interactions as described by Schunk and Nagy [2009]. 
Further information about this calculation is provided in 
section A2. Using IRI for the mean ion mass, we derive v in 
as shown in Figure 8b, which fits an exponential model with 
a scale height of 9.3 ± 0.7 km and a collision frequency of 


170 ± 10 s -1 at 119.9 km. We emphasize that the data 
points in Figure 8a are from observations only, whereas the 
data points in Figure 8b are derived in part from the IRI mass 
height-profile. 

3.1.2. Validation of in Situ Measurements 

[ 16 ] Since the ion-neutral collision frequency is propor- 
tional to neutral density, the re scale height H K is a good 
approximation for the neutral density scale height over small 
ranges of altitude. To estimate the height variation of H K , we 
performed least squares fits of the exponential model to 
successive pairs of re measurements between 117 km and 
131 km, the region over which re varies approximately expo- 
nentially. Results are shown in Figure 9 which includes the 
NRL-MSISE-00 model scale height (dashed line) calculated 
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Figure 6. Time series plots of total SI1 signal, with four- 
parameter Gaussian fits (dashed curves). The gray curve in 
each panel is the component of the E-cross-B drift velocity 
projected into the plane of the SII aperture. The labeled arrows 
show the equivalent phase shifts between the SII signal and the 
E-cross-B drift in the rammed payload frame, (a) An example 
of the profiles used for calibrating the E-cross-B direction and 
the measurement uncertainty, (b and c) Zoomed examples 
from the profiles used to calculate re. (d) An example profile 
for which equation (1) has no solution. 
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Ion ram deflection (36.234) 



Figure 7. Angular deflection of ion ram vector (v ' ± ) from 
the E-cross-B vector (v BxB ) as observed in the payload refer- 
ence frame. Below 117 km the measured angular deflection 
is anomalously large and we are unable to obtain estimates 
for re below this altitude. 


along the rocket trajectory for reference. Neutral atmospheric 
pressure scale heights (solid line) calculated from the ioni- 
zation gauge measurements, using MS1S for the neutral 
composition, are also shown. The scale-heights derived from 
the S1I ion ram deflection agree very well with ionization 
gauge (1G) neutral scale height measurements below 125 km, 
providing a convincing validation of the re measurements in 
this region. 

3.2. PFISR Method and Results 

[ 17 ] Typically, incoherent scatter spectra or autocorrela- 
tion functions are fit for electron and ion temperatures (as 
well as other parameters, like Doppler motion and electron 
density) under the assumptions of a composition profile and 


a collision frequency. In the lower E region (90-100 km), 
the plasma is collisional and it is reasonable to assume that 
electron and ion heating are negligible and that the plasma is 
in thermal equilibrium with the neutrals. In this collisional 
regime, the 1SR spectrum becomes Lorentzian in shape [e.g., 
Dougherty and Farley, 1963] and the correlation time 
(inverse of spectral width) is proportional to vT l , where v 
is the momentum transfer collision frequency and T is the 
temperature. Thus a spectral change is ambiguous in terms 
of a change in collision frequency and temperature. For the 
purposes of IS fitting, this ambiguity makes it impossible 
to fit for temperature and collision frequency at the lower 
E region altitudes. At slightly higher altitudes, it may be 
possible to fit for collision frequency and temperature simul- 
taneously, however one must be cautious of the breakdown 
in assumptions about thermal equilibrium as the ions are 
heated by frictional heating and the electrons are heated by 
wave-particle interactions [e.g., St.-Maurice et al., 1981; 
Hysell et al., 2009], 


(a) K=El c !v\ n (36.234 descent) 




Figure 8. Height profiles of (a) re, the ratio of ion cyclotron 
angular frequency to ion-neutral momentum transfer colli- 
sion frequency v in , and (b) the ion-neutral moment transfer 
collision frequency. The solid curves are least squares fits 
to exponential models. The long-dashed curves are the theo- 
retical results calculated according to the prescription in 
section A2. In Figure 8a, the data points are calculated from 
observations only, whereas in Figure 8b the data points are 
derived from the observed re using the 1RI model for the 
mean ion mass. 
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Figure 9. Height profile of the neutral density scale height 
derived from measurements of neutral temperature (solid 
line) and measurements of the k scale height (points). Also 
plotted is the MSISE-00 neutral scale height (dashed line) 
for rocket 36.234. 

[is] An alternative technique is based on measurements of 
the variation of ion velocity with height, using equation (1) 
to obtain k. This approach is similar to the analysis of the 
rocket data, except that in this case the full ion velocity drift 
vector perpendicular to the magnetic field is measured, and 
there is no independent measurement of electric field. To 
obtain n we 1) resolve line-of-sight velocity into a vector as 
a function of altitude; 2) determine the E x B drift from data 
at altitudes greater than 200 km; and 3) solve equation (1) 
for k as a function of altitude using the same approach as 
Sangalli et al. [2009], with and without the inclusion of the 
neutral winds measured by the rocket. Individual measure- 
ments are resolved at 5 -minute resolution over two 15 -min 
intervals, centered on the two TMA rocket launches: 
12:24:54-12:40:13 and 12:40:14-12:55:32. We estimate k 
from the mean of those three measurements, and estimate the 
uncertainties as the standard deviation over y/3. The data and 
results are shown in Figure 10. Shown are observations and 
results derived using winds from TMA rocket 41.064 laun- 
ched at 1231 UT (Figure 10, top) and TMA rocket 41.065 
launched at 1246 UT (Figure 10, bottom). Figure 10 (left) 
shows vector winds (smooth curves), smoothed over 8 km, 
and ion velocities (points) in geomagnetic coordinates, 
where black is East and gray is North. Figure 10 (right) 
shows n derived with (black) and without (gray) the neutral 
wind taken into account. Least squares fits (dashed lines) 
using equation (4) give z 0 = 115.9 ± 1.2 km (top) and 
z 0 = 119.9 ± 1.1 km (bottom), with scale heights of 
9.1 ± 1.0 and 7.3 ± 0.9 km, respectively. 

[ 19 ] Figure 11a compares the in situ measurements from 
rocket 36.234 (dots) with those from the PFISR calculated 
from averages over the interval 12:24:54 UT to 12:40: 13 UT 
(squares). Figure lib compares the in situ measurements of 
Sangalli et al. [2009] with those of PFISR averaged over the 
interval 12:40:14 UT to 12:55:32 UT. 

4. Discussion 

[ 20 ] The smoothness of the ion signal within a payload 
spin period suggests that there is very little structure in either 


the ion density or ion temperature over the course of a 
payload spin period, which corresponds to spatial scales less 
than 1 km. Large electric fields were observed on the upleg 
portion of 36.234, but the downleg had field strengths less 
than 5 mV/m, consistent with 7) — T n < 20 K \Saint-Maurice 
and Hanson, 1982], assuming that there are not significant 
wave-particle interactions. This validates one of the key 
assumptions that goes into the analysis of PFISR spectra in 
this altitude region, at least for this case study. 

[ 21 ] Often it must be assumed that neutral winds are neg- 
ligible when estimating collision frequencies from ISR data. 
In our case study, comparison of the PFISR n profiles that are 
calculated with and without neutral winds is consistent with 
previous work [e.g., Thayer, 2000] that shows ignoring 
neutral winds can introduce significant errors on estimates of 
k (and hence collision frequency) in the transition region. 

[ 22 ] Another assumption that goes into the analysis of low- 
altitude PFISR data is that the ionosphere is homogeneous 
and static over scales and times sampled by the radar beams. 
In situ observations can provide insight into the validity of 
this assumption for the event studied here. During the first 
interval (Figure 8), between altitudes of 1 18 km and 128 km, 
the in situ measurements show a 20% decrease in k relative to 
the empirical model. This difference is consistent with a local 
enhancement in neutral density relative to the model. This 
observation provides evidence for localized Joule heating 
below 125 km causing a localized vertical (and/or horizontal) 
structuring of neutral mass density over scales of a few km. 
The in situ observation ofz 0 reported here is 1.9 km higher 
than that reported by Sangalli et al. [2009] for the period 
approximately 16 min later, and 4 km higher than that 
obtained by the radar. In the case of the two in situ mea- 
surements, the difference cannot be explained by the reported 
measurement uncertainties, so it is due either to unaccounted- 
for systematic errors in one or both of the measurements, or it 
is a true difference. Due to the exponential falloff in neutral 
density with height, a difference of 1.9 km in z 0 represents a 
significant variation in neutral atmospheric structure. Also 
significant is the large difference in the atmospheric scale 
height obtained from these two rockets. Examination of 
Sangalli et al. [2009, Figure 8] reveals a relatively constant 
scale height of about 5 km in the range 108 km to 128 km, 
about half the value reported here. The two rocket flights were 
16 min apart and their lateral separation was approximately 
1 1 km at a height of 120 km on their descents (see Figure 4). 
There is an excellent agreement between the radar and in situ k 
measurements for the second interval (rocket 21.138) over a 
limited height range of 108-1 18 km (Figure 1 lb). There are 
significant differences outside this altitude range, and the in 
situ estimate for z 0 during the second interval (119. 9 ± 1.1 km) 
does not agree with the in situ estimate from Sangalli et al. 
[2009] (118.0 ± 0.3 km) at the 1 -sigma uncertainty level. 
Piecing these observations together paints a picture of an 
atmospheric neutral density that is structured at scales of 1- 
1 0 km, vertically and/or horizontally, and that is also changing 
in time. This structure might explain some of the differences 
between the radar and in situ observations, since the radar 
measurements are averaged over larger distances (on the order 
of several km horizontally) and longer times (averages of three 
five-minute integrations). Notwithstanding this qualification, 
incoherent scatter radars can provide accurate measurements 
of IT coupling. We note the excellent agreement between 
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2007 - 01-19 12 : 24:54 - 12 : 40:13 



2007 - 01-19 12 : 40:14 - 12 : 55:32 



Figure 10. Determination of re from PFISR observations. Shown are observations and results derived 
from two 15-minute intervals centered on the TMA rocket launches: (top) TMA rocket 41.064 launched 
at 1231 UT and (bottom) TMA rocket 41.065 launched at 1246 UT. (left) Vector velocities of winds 
(smooth curves) and ions (segmented curves), where black is East and gray is North, (right) re derived with 
(black) and without (gray) the neutral wind taken into account. 


re scale height measurements made in situ (9.3 ± 0.7 km) and 
with the radar (9. 1 ± 1 .0 km) during the first interval. 

[ 23 ] Above 118 km the technique for determining re from 
a spinning spacecraft is robust since few assumptions go into 
the measurement. We have seen with rocket 36.234 that 
accurate estimates of the collision frequency based on 
deviations of the ion drift angle relative to the E-cross-B drift 


direction are only possible with precise calibration of the 
relative angle between the measurements to an accuracy of 
0.25°. This calibration takes into account misalignments of 
instrumentation and other sources of error. The results are 
insensitive to variations in spacecraft charging, which can 
affect measurements based on the magnitude of the ion 
velocity if not accurately accounted for. Moreover, there is 
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K=fl c /v j n : In situ and PFISR 



K=Cl c /v m : In situ and PFISR 



Figure 11. (a) Comparison of in situ (circles) and PFISR 
(squares) measurements of re obtained during the flights of 
the first pair of rockets 36.234 and 41.064. (b) Comparison 
of in situ (solid line) and PFISR (squares) for the period cov- 
ering the second pair of rockets 21.138 and 41.065. The 
solid curve in Figure lib, from rocket 21.138, is the data 
of Sangalli et al. [2009]. 

no assumption about the form of the ion velocity distribution 
function. In principle the in situ technique does not require 
calibration against electric field measurements when the 
angular position of the SI1 and electric field booms is known 
with sufficient accuracy and precision. Future experiments 
could be conducted using precision boom angle measure- 
ments to eliminate the need to calibrate the ion measurement 
in the region above 140 km. Variability in ion density and/or 
ion temperature within a payload spin period would, how- 
ever, reduce the effectiveness of the technique. 

[ 24 ] Below 118 km the deflection of the ion ram vector on 
rocket 36.234 is significantly larger than what is expected 
from the observed neutral wind and electric field observa- 
tions. An explanation may lie in the fact that the neutral wind 
measurements were obtained from the upleg of rocket 
4 1 .064, whereas the ion and electric field measurements were 
obtained on the downleg of rocket 36.234. Any variations in 
winds between the two times and different latitudes could 
affect the re estimates particularly at lower altitudes, where 
the ion drift is largely controlled by the neutral atmosphere. 

[ 25 ] Watcmabe et al. [1991] presented observations from 
three rockets of the ratio of v based on measurements of 


v,X; Ej_, and B, with the assumption that u± = 0, and found 
that above 125 km under quiet auroral conditions the ions 
drift according to 

Vi± =E ± x b/B (5) 

to a good approximation. Our in situ observations, from the 
downleg of rocket 36.234, were obtained from the north side 
of a weak auroral arc, and the weak electric fields on the 
downleg (<5 mV/m) are consistent with a quiet auroral 
environment. Deviations from the E-cross-B drift approxi- 
mation equation (5) are evident in Figure 7 at altitudes up to 
140 km, and show precisely at what altitude the approxima- 
tion can begin to break down under quiet auroral and geo- 
magnetic conditions. Watanabe et al. [1991] used the 
assumption u± = 0 and the measured scale height H and ratio 
vjili above 115 km to estimate the neutral wind altitude 
profile below 115 km assuming an exponential model for the 
collision frequency. For the present case, the MSIS and IRI 
models show that the collision frequency does not vary 
exactly exponentially over scales greater than a few km in 
this atmospheric transition region, but the observations from 
rocket 36.234 reveal that an exponential model for collision 
frequency appears to work well between 118 km and 130 km 
in this case. But as we do not have full ion velocity vector 
information from 36.234, we are unfortunately unable to 
estimate the neutral winds using the method of Watanabe 
et al. [1991], independently of the TMA measurements. 

[ 26 ] If the results from rocket 36.234 are taken as repre- 
sentative of quiet auroral conditions, we see that qualitatively 
the MSIS and IRI models do well in representing collision 
frequencies, but also that there are disparities between the 
empirical models and observations over scales of 1-10 km. 
Figures 8 and 9 show that, to within measurement uncer- 
tainties, the re scale height is approximately constant through 
the altitude region 118-130 km. Given that the neutral 
atmosphere is not isothermal in altitude, there is evidence that 
details of either or both of the average neutral particle and ion 
mass height profiles are not represented accurately at scales 
of 1-10 km by the MSIS and IRI models. These points serve 
as a reminder that the empirical models should be used with 
caution when applied to problems over scales of 1-10 km, 
even under quiet auroral conditions. 

5. Conclusions 

[ 27 ] 1 . This study presents a treatment and analysis of in situ 
(NASA rockets 36.234 and 41 .064) and incoherent scatter radar 
(PFISR) measurements to investigate ion-neutral momentum 
coupling of the high-latitude auroral ionosphere-thermosphere 
system. 

[ 28 ] 2. We have introduced a robust technique for esti- 
mating ion demagnetization in the E region ionosphere from 
an ion analyzer and an electric field antenna on a spinning 
spacecraft (Figures 6 and 7). 

[ 29 ] 3. Ion velocity and electric field measurements were 
combined with neutral wind measurements to estimate the 
height profile of re = £Ljv in between 118 km and 140 km 
(Figure 8) in a relatively undisturbed auroral E region. 

[ 30 ] 4. During the period 12:24:54 UT to 12:40:13 UT, the 
height of transition between the collisional and collisionless 
regimes (for which re = Q,/re,„ = 1) was measured in situ to be 
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z 0 = 119.9 ± 0.5 km (Figure 8), and by PFISR to be 
z 0 = 115.9 ± 1.2 km (Figure 10). The scale height mea- 
surements were in excellent agreement: 9.3 ± 0.7 km (in situ) 
and 9.1 ± 1.0 km (PFISR). Comparison of PFISR and rocket 
results for the period 12:40:14 UT to 12:55:32 UT show 
excellent agreement in the re measurements over a limited 
height range of 108-118 km, although the estimates for z 0 
do not agree to within measurement uncertainties. 

[ 31 ] 5. Measurements of the re scale height agree well with 
neutral atmosphere scale height measurements from an ion- 
ization gauge (IG) on the same rocket between 118 km and 
125 km (Figure 9). To within measurement uncertainties, 
the re scale height is relatively constant between 118 and 
130 km, significantly below the IG measurements at alti- 
tudes above 125 km. 

[ 32 ] 6. Variability in re for a single rocket, between rockets 
flown 16 min apart, and between in situ and ground-based 
observations, is evidence for substantial structuring of the 
atmospheric neutral density vertically and/or horizontally. 
This structure is suggestive of structured frictional heating at 
scales of 1-10 km. Such variations in space and time may 
account for the underestimate in z 0 by PFISR during the 
interval 12:24:54 UT to 12:40:13 UT. 

[ 33 ] 7. Qualitatively, MSIS and IRI models perform well 
in representing collision frequencies under relatively quiet 
auroral conditions, but there are disparities between the 
empirical models and observations over scales of 1-10 km 
(Figures 8a and 8b). 

[ 34 ] 8. There were no indications of significant structure at 
scales below 1 km in either ion density or ion temperature 
between altitudes of 118 km and 135 km (Figure 6). 

Appendix A 

Al. Uncertainty in n (in Situ Technique) 

[ 35 ] We estimate the uncertainty in re at a given height 
from 

A« = -AR (Al) 

where 6 is the deflection of the ram velocity vector from the 
direction it would have if the plasma were collisionless. The 
uncertainty in this angle, Ad, has two parts, 

A0 = \J (A8 m ) 2 + (A# e ) 2 , (A2) 

where the measured uncertainty is the error in the observed 
ram direction and the expected uncertainty is the error in the 
expected angle, which is derived from the field and neutral 
wind uncertainties using 



where the partial derivatives can be derived from the angle 
between v[ and vkxB (Figure 5). Uncertainties in each com- 
ponent of the neutral velocity vector are on the order of 
1-5 m s _I over the course of the ten minute observation 
period. Uncertainties in the electric field components are of 
order 0.5 mV/m. 


A2. Collision Frequency Calculation 

[36] The momentum transfer collision frequency height 
profile is calculated using the NRL-MS1SE-00 model mass 
densities for O, N 2 , and 0 2 , and the dominant IRI ions. The 
A p and F l01 averages for the period leading up to the flight 
are included in the model. The collision frequencies are 
obtained for all pairs of ion-neutral interactions from the 
polarization attraction collision frequency from equation 
(4.88) of Schunk and Nagy [2009]. Total collision frequen- 
cies between an ion and all other neutral species are obtained 
from a sum over all neutrals, and the average momentum 
transfer collision frequency is a density-weighted average of 
the total collision frequencies for all ions. Relative ion con- 
centrations for NO + , 0 2 , and 0 + are obtained from IRI 2007. 
All other IRI species have negligible concentrations at these 
altitudes. Uncertainties in the calculated collision frequen- 
cies can arise from uncertainties in ion or neutral concen- 
tration, and ion or neutral temperature. Neutral temperature 
is obtained from NRL-MSISE-00, and the ion temperature is 
assumed to be equal to the neutral temperature, which is 
justified given that the electric fields on the downleg were 
less than 5 mV m _1 . 
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